ON THE EFFECTS OF BOHM'S POTENTIAL ON A MACROSCOPIC 
SYSTEM OF SELF-INTERACTING PARTICLES 



OLIVER TSE* 

Abstract. Wc consider an instationary macroscopic system of self-interacting particles with an 
additional potential, the so called Bohm's potential. We study the existence of non-negative global 
solutions to the (4-th order) system of equations and allude the differences to results obtained for 
classical models. The problem is considered on a bounded domain up to three space dimension, 
subject to initial and Dirichlet boundary condition for the self-interacting potential and Neumann 
boundary condition for the particle density. Moreover, the initial datum is only assumed to be 
non-negative and to satisfy a weak integrability condition. 
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1. Introduction. Consider a macroscopic system of self-interacting particles 
with Bohm's potential, which describe the evolution of the normalized density n > 0, 

(la) n t = div(nVF) in O, nd v F = on V, 

with the quasi fermi-level F given by 

(lb) F = -e 2 ^^ + logn-cr$ in 17, d v \fn = on T, 



and the potential $ due to self-interaction in a particle system, 
(lc) -A$ = ninft, $ = on T, 

supplemented with the initial condition n(0, ■) — uq > 0, where v is the outer normal 
to the convex, bounded domain SI C R d , d < 3 with Lipschitz boundary T, e > is the 
scaled Planck constant and a £ R is the mass of the system of self- interacting particles. 
In this case, positive mass a > would indicate the presence of self-attraction, while 
negative mass a < self-repulsion. Observe that the presence of the Bohm potential 
yfn Ay/n, which is a non-local second-order term, leads to a fourth-order evolution 
equation for the normalized density n, given by 

(2a) m = div ^nV ( - e 2 + log n - J mil, 

with the natural boundary conditions 

(2b) nd v { - e 2 ^^ + \ gn - er$) =0 on T, and d v \fn = on T. 



Clearly, the boundary conditions for n imply J„ n dx — J n hq dx. Therefore, it is 
sufficient to ensure that no dx = 1 for n to be kept normalized. 
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Such kind of systems occur typically in the theory of transport equations for 
semiconductors as a macroscopic limit (quantum drift-diffusion equations) of its mi- 
croscopic counterpart (Wigner-Boltzman-Poisson system) (c.f. [13, 16] and references 
therein). In this setting, n describes the electron density and therefore o < 0, since 
electrons are negatively charged. In the case of no interactions, F = —^fn A-^/n, 
also called the Derrida-Lebowitz-Speer-Spohn (DLSS) equations, the authors in [10] 
and [15] successfully proved the existence of solutions and additionally showed the 
rate of convergence of the solutions to the unique stationary solution. Many of the 
ideas in this paper follow as a result of these two papers. 

By passing to the limit e —¥ 0, we formally recover either the classical drift- 
diffusion equations (a < 0) or the standard model for a system of self-gravitating par- 
ticles (a > 0) [3]. This system is well known to have global solutions for a G (— oo, 87r) 
and blowup phenomena for large mass a > 8ir in d = 2. Another model of the limit 
system comes from statistical mechanics for vortex points as the mean field limit 
of the canonical Gibbs measure associated to a iV-vortex system in a bounded two- 
dimensional domain [5,6]. Here, the authors established the concentration phenomena 
for the weak limits of the Gibbs measures, when N — > oo, to the solution of the limit 
system (e = 0) in the case a G (— oo,87r). Similar systems that manifest blowup 
phenomena arises in the theory of combustion [2] and chemotaxis equations [24]. 

The purpose of this paper, among other things, is to show that the presence of 
the Bohm potential (e > 0) leads to a regularization of the limit problem (e = 0). 
More precisely, it is shown that (1) possesses at least one global weak solution for any 
a G R, i.e., the presence of the Bohm potential prevents any blowup phenomena. A 
short discussion on stationary solutions to (1) for a G R can be found in [19]. There it 
is also shown that stationary solutions with e > 0, d = 2 and a G (—oo, 8tt), converge 
weakly to stationary solutions of the limit problem (e = 0) . 

For notational convenience, we set 

V := {u G H 2 (n) | d v u = on V} , V := H 2 (n) n fl£(fi), 
and denote the space of probability measures with finite entropy £ by 

V := ju G L^fi) | u > 0, J udx = 1, £(u) < +oo| , 

where £ is simply the (negative) physical entropy given by 

£(u) := / (ti(logit-l) + l)da:. 
Jn 

The main purpose of this paper is to show the following result. 
Theorem 1. For any finite T > and n G V, there exists at least one solution 
(n, F, $) to system (1), with 

neW 1,1 (0,T;V*), y/n G L 2 (0, T; V), 
FVneL 2 (0,T;I 2 (O)) and $ G L\0, T; Vq), 

and additionally n>0 a.e. in (0,T) x O. 

The proof of Theorem 1 relies on the fact that the physical entropy £ provides a 
controlled growth estimate for (2). Indeed, by formally multiplying (2a) with log(n), 
integrating over f2 and integrating by parts, we obtain 

d f 9 ,2 

— £(n)+c / |V 2 \/n| da; < ci, 

di Jn 
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for some constant c , C\ > 0, independent of n. By using interpolation inequalities on 
derivatives [1], i.e., there exist constants 5 > and c(<5) > such that 

\\Vu\\ 2 2 <S\\V 2 u\\t + c(S)\\u\\l 

along with the logarithmic-Sobolev inequality [23] 

/ u 2 log(-^)dx <c L \\Vu\\ 2 2l 

where cl > only depends on and d, and the fact that ||n||i = 1, we further obtain 

^f(n) + c 2 £(n) < c 3 , 

for some constants 02,03 > 0. An application of Grownwall's lemma on this estimate 
provides the global boundedness in time t of £{n), and consequently shows (formally) 
the absence of a blowup phenomena in the space V . 

2. Preliminary results. In this section we recall several standard results re- 
garding unique solvability and regularity of solutions for elliptic equations. We be- 
gin by recalling a well known interpolation theorem for Sobolev spaces, namely the 
Gagliardo-Nirenberg-Sobolev inequalities [26]. 

Proposition 2. Let O C M d be a bounded Lipschitz domain, m e N and 1 < p, 
q, r < oo. Then there exists a constant c > 0, such that 

\\D a u\\ r < c WuWljuW 1 - 6 for any u G W™(ft) H L«(Sl), 

provided that < \a\ < m — 1, = \a\/m and \a\ — d/r = 8(m — d/p) — (1 — ff)d/q. 
If m— \a\ — d/p ^ No, then the values \a\/m < < 1 are allowed. 

We also recall a regularity result for linear elliptic problems on convex, bounded 
domains due to [12]. 

PROPOSITION 3. Let £1 C R d be a convex, bounded domain and f E L 2 (£l). Then 
the homogeneous Dirichlet problem 

— Au = f in ft, u = onT, 

possesses a unique weak solution u € Vq, which satisfies the estimate 

N| 2 ,2 < c||Au|| 2 , 

for some constant c > 0, depending only on the diameter of Q,. 

The next result we want to recall is a regularity result for linear elliptic equations 
with homogeneous Neumann boundary conditions and source terms from the Orlicz 
space LlogL(f2) (c.f. [4, 22]). For the sake of completeness, we include a simple proof 
for this result in Appendix A. 

PROPOSITION 4. Let fl cM. d be a bounded Lipschitz domain, f e LlogL(fi) and 
a e L°°(Q) with a > a > a.e. in ft. Then the elliptic equation 

(3) — div(aVu) = / infl, d u u = Q onT, 

has a unique weak solution u € Wp, with (1 + M) 1 / 2 € H 1 ^), where 

W := j« e W h ^{Sl) | Judx = f3^ with \\u\\ Wfi :=\\Vu\\_i_, 
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and ft £ K is some given constant. 

Another result we will need is a result regarding nonlinear elliptic equations with 
natural gradient growth [19] (c.f. [8, 9]). 

Proposition 5. Let ft C R d be a bounded Lipschitz domain and f £ L P (VL), 
p > d/2. Then the elliptic equation 

-y (^Ay+i|Vy| 2 ^ +y = f m ft, 5„y = on T, 

has a weak solution y £ if 1 (ft) n L°°(Q). Furthermore, e«/ 2 e ff^ft) n L°°(ft). 

3. Existence of solutions. The strategy for proving the existence of solutions 
involves first showing existence of solutions for a time-discrete problem with time step 
t > with the help of Leray-Schaudcr's fixed point theorem. By establishing uniform 
bounds on the sequence of solutions {n^} with respect to t then allows to extract a 
subsequence, which converges to a solution of (1) when passing to the limit t — > 0. 

3.1. Existence of a time-discrete solution. Let r > an arbitrary but fixed 
time step and w £ V be a given function. The task at hand is to find a weak solution 
(F, s/n, $) £ Wp x V x Vo, for some /? € R, to the semi-discrete equation 

(4a) -div(nVF) = (to - n)/r in f2, n^F = on T, 

(4b) -e 2 ^^+logn = (T$ + F in 0, d„y/n = on T, 

(4c) -A$ = n in 0, $ = on T, 

where F, n and $ satisfy (la), (lb) and (lc) respectively. 

We mention some of the problems encountered in solving this problem. Observe 
that, in order to solve (4a) for F, we have to make sure that n £ L°°(f2) is uniformly 
bounded below away from zero. Otherwise, we will have to rely on degenerate elliptic 
estimates, which will require some kind of regularity on both n and n^ 1 . 

To circumvent both these problems simultaneously, we introduce a transformation 
for n, given by n = e y /\\e y \\i. Clearly, if y £ L°°(Cl), then n satisfies the assumptions 
in Proposition 4. Moreover, the normality of n is satisfied trivially. Inserting this into 
(4b), leads to a nonlinear elliptic equation with natural gradient growth, 

(5) -^Ay + ^\Vy\ 2 ^ +y = a<i> + F + \og\\ e y\\ 1 inn, d v y = on T. 

We will then use (5) to construct an auxiliary problem, similar to (4), and apply the 
Leray-Schauder theorem on this problem. Under the appropriate regularity, we may 
then recover solutions of the original problem (4) . 

Lemma 6. Let d < 3 and w £ V . Then there exists at least one weak solution 

(F, y/n, <&) £ Wp x Vx Vo, 

to (4) with /3 = log He* ||i and n£V. 

Proof. Let w £ V . For arbitrarily given v £ C(f2) and A £ [0, 1], wc consider the 
auxiliary problem to find (y,F,&) £ V x Wp x Vo, with (3 = -A|Q| log ||e"||i: 

(6a) -div((e7||eli)VF) = \(w - (e"/||e"||i))/T in ft, d u F = on F, 
(6b) -^Ay + i|Vy| 2 ^) + ?y = a$ + F + Alog ||e"||i in ft, d v y - on F, 
(6c) -A* = A(e^/||e^ || x ) in ft, <f> = on T, 
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As mentioned above, we used the variable transformation n = e y /\\e y \\i. Note also 
that for the case A = 1, system (6) is equivalent to the initial system (1). 

We begin by showing that this system yields a weak solution y G V for any 
v G C(Q) and A G [0,1]. Indeed, since e"/||e"||i G C(fl) and is uniformly bounded 
from below by a positive constant, we obtain a unique solution F G Wp for (6a) due 
to Proposition 4. From Proposition 3, we obtain a unique solution $ G Vq for (6c). 
Now, by inserting (F, $) into (6b), we obtain a solution y G i? 1 (fi) n L°°(Q), due to 
Proposition 5. Note also that n G L°°(fl) with n > e iafi/ / 1 1 e i/ 1 1 x =: R > a.e. in Q, and 
J Q ndx = 1. Furthermore, by rewriting (6b) in terms of p := y/n, we have 

-e 2 Ap = p(a$ + F- log p) e i 2 (^), 

which implies Ap G L 2 (f2), since p G L°°(Q). Due to the convexity of fi and the 
homogeneous Neumann boundary condition for p, we have the second order inequality 
||V 2 p|| 2 < ||Ap|| 2 , where V 2 p denotes the Hessian of p [12]. Therefore, p e # 2 (f2) 
and consequently j/ G V, since 

(7) %y = 2d l [ p - 1 d 3 p) = 2(pd ijP - d i pd jP )/p 2 e l 2 (0), 

Due to the boundary condition for y, we further have d v p = a.e. on T, i.e., p G V. 

Now consider the operator H:C(f2) x [0, 1] — > C(S1); (u, A) y. This operator is 
continuous and compact due to the continuity of the solution operators corresponding 
to (6a)-(6c) respectively and the compact embedding H 2 (£l) C(f2). It is also easy 
to see that H(u, 0) =0 for all v G C(Q). We see this by simply testing the variational 
formulation of (6b) with tp — sign(y)(el y — 1), which yields 



/ \Ve^' 2 \ 2 dx + ^- [ \Vy\ 2 dx + f \y\ 2 dx < 0, 
Jn 4 J n J n 



where we used the fact that s(e s — 1) > s 2 for s > 0. Therefore, y = a.e. in Q and 
consequently y = in f2, since y is continuous. 

Let (y, A) G C(£l) x [0, 1] such that H(y, A) = y. We now show that y is uniformly 
bounded (in A) by some constant M > in C(fi). Observe that for -^n G V, the 
auxiliary system (6) is equivalent to the equations 

(8a) -(n-w) = div fnV ^ , + ] () „ „ _ 



(8b) -A$ = An, 



where we also used the fact that log ||e y ||i is constant. Since (j)(s) = s(logs — 1) + 1 
is convex, <p(si) — 0(s 2 ) < (p'(si)(si — s 2 ) for all si,s 2 > 0. Therefore, 

A A f A f 

— {£{n)—£{w)) = — I (4>(n) — 4>{w)) dx < — / (n — w)logndx 
t T Jn T Jn 



2 A v / n 



nV — e — — + log n — cr$ • V log n dx 



= e 2 [ nV -Vlognda;- / n|V log n\ 2 dx + a [ W^-Wndx 

Jn V V n J Jn Jn 



= e Ji - J 2 + crJ 3 . 
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For the following computations we set p — y/n. For J\, we have the following 
Ji = 2^ (v(pAp) -2ApVp) - ^ydx 



= -2 f |Ap| 2 + Ap^dx<-2 / 
Jn P Jn 





V 2 p 


\^ 


Vp 


2" 




p 


p 


p 





dx, 



where, in the last inequality, we used the fact that || V 2 p|| 2 < ||Ap|| 2 , due to the results 
in [12]. As for J 2 , we have the Fisher information 



■h 



n|Vlogn| 2 dx = 4 



/' 2 

Jn 



Vp 



dx. 



Since $ satisfies (8b), we have 



/ (A$)ndx + / n<9„$ds = A / n 2 dx + / nd v 
Jn Jr Jn Jr 



$ds. 



Due to the regularity of $ and the boundary condition $ = on T, we have d v Q < 0. 
Therefore, J3 < J Q n 2 dx for any A € [0, 1]. Altogether, we have 



-(£(n)-£(w))<-2 f p 
T Jn 

From Appendix B, we obtain the estimate 

[ n 2 dx < Si [ 
Jn Jn 



e 2 


V 2 p 


2 Ap 
+ e 2 — 


Vp 


2 

+ 2 


Vp 


2 

a 

-~2 n 




P 


P 


P 




P 



dx. 



Vp 


4 

dx + 5 2 




Vp 


p 




in 


P 



dx + 0(61,82, a), 



for any 61,62 > and some a e (0, 1) for d = 2 or a e (1/3, 1) for d = 3. By choosing 
62 = 4/l°i and rescaling <5i = 2e 2 8/\a\ for some (5 > 0, we further obtain 



(9) 



(£{n)-£(w))+2e 2 J p 2 
Jn 





V 2 p 




Vp 


2 


Vp 


4" 






P 


p 




p 





dx < c{6). 



In Appendix C, we show that the second term on the left is bounded from below by a 
multiple of ||V 2 p||2, i.e., we have for some 7 > (explicitly given below) the estimate 



(10) ~(£(n) -£{w)) +2 7 e 2 / | V 2 p| 2 dx < c(5) for any A e [0, 1], 

T Jn 

which gives the uniform bounds 



(11) 



£{n) < £{w) +c(6)t and ||V 2 p||^ < 



1 



2 7 e 2 r 



(£(w) + c(8)t) , 



thus implying the boundedness of V 2 p in L 2 (f2) independent of A e [0, 1], and there- 
fore the uniform boundedness of p e V. As a consequence, we obtain an upper bound 
for y, i.e. there exists some constant M u > 0, such that y < M u a.e. in il. 

To show that y is uniformly bounded from below, we have to show a uniform 
lower bound for p away from zero. This is a result of the Harnack inequality [20, 25]. 



On the effects of Bohm's potential on a macroscopic system of self-interacting particles 



7 



Firstly, note that logp — er^ — F G L p (fl), p > d/2. Then by Harnack's inequality, a 
weak solution p G Hl oc (Q) of the equation 

-e 2 Ap + pp = in 0, with p G L p (fl), p > d/2, 

in a cube Q = Q(3r) C O with < p < M u in Q, satisfies 

Q(r) rV - Q(r) FV 

for some constant Ch = c/j(A) > independent of p. 

Now it is an easy exercise to verify that p(x) for some x G f2 would lead to p = 
in fi, which clearly contradicts j n p 2 dx = 1. Therefore, p > p > a.e. in il uniformly 
(in A). Hence, there exists some constant Mi > such that y > —Mi a.e. in Q,. 
Choosing M = max{Mj,M u } gives the estimate ||y||oo < M. Furthermore, we see 
from (7) that V 2 y G L 2 (il) is uniformly bounded. This implies that y is uniformly 
bounded in H 2 (Q), and consequently in C(Q) with the same constant M, due to the 
Sobolev embedding H 2 (fl) C(S1). We finally conclude the proof by applying the 
Leray-Schauder fixed point theorem [11]. □ 

Remark 7. The value |3€t may be thought of as a Lagrange multiplier for the 
constraint on the density J^ndx = 1. Indeed, the solution ^fri G V of (4b) may be 
characterized as the unique minimizer of the functional 

F{n) :=e 2 [ \\7 ^fh~\ 2 dx + E{n) - [ n(a<P + F) dx, 
Jn Jn 

on the space {n G "P | \fn G V} (c.f. [7]). 

3.2. A priori estimates. Let T > be a fixed arbitrary terminal time. For 
every r > we define the step function n^: [0, T) — > L X (Q) recursively as follows. Let 
n(0) = no and for given k G N, let (Ft, ^/n~k, $&) G Wp x V x Vo be the non-negative 
solution of (4) with w = nt-i- Now define 

F^{t) :=F k , n {T \t) := n k , $ (r) (t) := $ fc for (k - l)r < t < kr. 

Then (fM,nM,$( T )) satisfies 

(12a) (n (T) -Crn (r) )/^ = div(n (r) VF (T) ) in Q, 

(12b) fM - _ £ 2^^1 + lQgn (r) _ ; n Q 

(12c) -A$ (r) = n (r) in 0, 

together with their respective boundary conditions. Here, ( T denotes the shift operator 
(C T n {T) ){t) = n^(t - t) for t G [t,T). As a consequence of Lemma 6, we obtain the 
following uniform estimate for the sequence of step functions {n^}. 

Lemma 8. There exists a S > such that for any 5 G [0,<5oL the sequence of 
step functions {n^} satisfies 

i||n (r) - Crn (T) ||Li+^T,T;V*) + II ™ (T) II H + * (0,T;H^ (Q)) < C 

for some constant c > independent of r. 
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Proof. For notational convenience, we set as before p^ = V n( T ). We begin by 
establishing a uniform (in r) bound for the translations. To do so, we multiply (12a) 
with an arbitrary function <p e V, integrate over and integrating by parts to obtain 

(13) - f (n^ -( T n^)tpdx = [ F< T >(Vn< T > ■ V tp + Atp) dx 
T Jn Jn 

< (\\F^Vn^\\,+\\F^n^\\ 2 )M v , 
where we used the embedding H 2 (fl) =-> W 1,4 (Q). This implies the estimate 

- Crn (T) |li,tt( T ,r ; v.) ^ c i ^ ll^Vn^H^ + ||FW(t)nW(t)||J +4 dt, 
for some constant ci > 0. Since there are constants a > 0, such that 

\\FVn\\ 1 + 5 < c 2 \\Fp\\ 2 2 + c 3 \\\7 P \\f> W and ||Fn|| 2 < c 4 ||i^ + c 5 |M|^), 

3 

with 77(5) = (1 + 5)/(l — S) > 1, we further obtain 

(14) lll^-C.^H^^^^ < c 6 ^ ll^ll^ + IIVp^ll^ + H^V^Illdi. 

Therefore, the required estimate follows from the uniform boundedness of the right 
hand side of the above equation in r. 

We now make several observations from the results obtained in Section 3.1. As a 
direct consequence of the conservation of mass, we obtain the uniform bound 

\\p (T) \\L~(o.T;L 2 (n)) = II™ (t) IIl~(o,t : li(q)) = / n dx = l. 

Jn 

Furthermore, from (10), we deduce the uniform bound 

(15) ||V 2 p M || i2(0 ,T;L^)) <C , 

with a constant cq > depending only on d, fl, T, and no, and consequently, the 
uniform boundedness of the sequence {p^} C L 2 (0,T;V). In order to establish 
uniform boundedness of the first two terms in (14), we make use of Proposition 2. For 
the first term, we have 

(16) £ \\p\\ 2 ^ dt < c 7 \\ P \c-$$ m £ \\p\\% {5) dt, 

with 8 = d/A e (0, 1) for d < 3, which leads to uniform boundedness if we choose 
T}(S) < A/d. This is equivalent to the requirement 6 < (4 — d)/ (4 + d) =: 61. Similarly, 
we have for the second term 

(17) £ || VHlf {5) dt < c 8 \\p\\lt-$$ m £ \\p\\ 2 2 9 f 5) ^, 

with 6 = (4 + d)/8 G (0, 1) for d < 3. Here, we may choose T](5) < 8/(4 + d), which 
is equivalent to choosing 6 < (4 — d)/(12 + d) =: S 2 . Therefore, we may choose 
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So = min{(5i, 5 2 }- The uniform boundedness of the last term in (14) may be seen as 
follows. From (12b), we obtain the following estimate 

\\Fp\\ 2 < e 2 ||V 2 p|| 2 + 2||plogp|| 2 + M||p$|| 2 < e 2 ||p|| 2 , 2 + (2 + c 9 \<j\)\\p\\l 

where we used the fact that slog(s) < s 2 for s > 0, the a priori estimate provided 
by Proposition 3, and the Sobolev embedding i^ 2 (SI) =-> L°°(£l). Now, from Proposi- 
tion 2, we deduce the existence of constants c w ,Cn > 0, such that 

(is) ll^/^lli < c 10 ^ ll/^lll. c i=7o^;^-c^>> H^ll^^ cl*, 

with 9 = d/8 £ (0, 1), which shows the uniform bound for {F^ p^}. Therefore, by 
collecting the estimates obtained above, we conclude the first assertion. 

To show that {n^} is uniformly bounded in L 1+s (0, T; H 2 (fl)), we simply use 
the fact that V 2 n = 2(|Vp| 2 + /oV 2 p) and Proposition 2 again, to obtain 

f \W 2 n\\l +s dt < c 12 f || vHlf +5) + IIpII 2 ^ + II v 2 HI^ at. 

Jo Jo 

Since (1 + 5) < r)(6), we have the uniform boundedness of the right hand side due 
to the estimates (15), (16) and (17). Similarly, we can show the uniform bound for 
{n^} in L 1+s (0,T;L 2 (fl)), which leads to the estimate asserted for {n^}. □ 

3.3. Passing to the limit r — > 0. We begin this section by recalling a general- 
ized version of Aubin's lemma found in [21]. 

PROPOSITION 9. Let X,B,Y be Banach spaces and X ^ B ^ Y with compact 
embedding X B, and F C L q (0,T]B) n £*((), T;X) ; q e (l,oo] ; be bounded with 

\\u — (ru\\Li( Tt T:Y) ^0 as t — > 0, uniformly for iieF, 

then F is relatively compact in L p (0, T; B) for all p < q. 

A simple consequence of Proposition 9 is the following result. 

Lemma 10. There exists a nonnegative function n <G M /1 ' 1 (0, T; V*) such that the 
following convergences hold true for some subsequence of {n^}: 

(n< T > - Crn (T) ) /t -» d t n in L 1+s (0, T; V*), 

n (r) ^n and Vn' T ' -> Vn a.e. in (0,T) x O, 

for t — > ; where — denotes the weak convergences in their respective spaces. 

Proof. Since the Lebesgue-Bochncr spaces L 1+<5 (0, T; V*) and L 1+5 (0, T; H 2 {n)), 
with S > are reflexive, they are weakly sequentially compact. Due to Lemma 8, 
we obtain a subsequence of {n^} (not relabeled) and some n € L 1+s (0, T; H 2 (fl)) 
such that the second convergence holds true. Furthermore, Proposition 9 provides yet 
another subsequence of {n^} such that 

n (T) ->n inL 1+4 (0,T;W ll4 (fi)) ) 

due to the compact embedding H 2 (Q) VF 1,4 (f2). In particular, we may extract 
a subsequence such that — > n and Vn' T ' — > Vn almost everywhere. Moreover, 
since n/ 1 "' > uniformly almost everywhere for all r, the limit function also satisfies 
n > almost everywhere. 
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In a similar way, we obtain some \ G L 1+s (0, T; V*) such that 

(ii (t) - Cr^ (T) ) /t^x in (0, T; V*) 

for some subsequence. In order to identify this limit with the time derivative of n, we 
multiply n with arbitrary functions v £ V and tp € Cq°(0, T), integrate over time and 
space, and integrating by parts w.r.t. time, we obtain 

[ T {n,v)vvtdt= lim /^V^) ^'^^ dt 
Jo r ^°Jo r 

= lim / ( ^ ,v) v <pdt = { X (t),v)v<p(t)dt. 

t ^°Jt t Jo 

Since the set {<pv \ v e V, ip e Cg°(0,T)} is dense in L"(0,T;V), q = (1 + S)/S, we 
have by definition of the generalized time derivative that dtn = □ 

Lemma 11. There exists a nonnegative function p = \fn € L 2 (0,T; V) such that 
the following convergences hold true for some subsequence of {p^}: 

p {T) ^ p mL^{Q,T;W^\VL)), 
p (T) -> p and Vp (r) -> Vp a.e. in (0,T) x ft, 

for t — y 7 with the same S > as in Lemma 10. 

Proof. Due to the boundedness of {p( T >} in L 2 (0,T;V), we obtain a weakly 
convergent subsequence (not relabeled) and some p such that ^ p in L 2 (0, T; V). 
Since — > n almost everywhere, so does p^ = V n( T ) -4- ^/ra almost everywhere, 
which allows us to identify p = ^fn > almost everywhere. Furthermore, we have 
that Vp( T -* —> Vp almost everywhere. Indeed, this follows from 

2p (r) Vp (r) = Vji (t) -^Vn= 2pVp a.e. in (0,T) x ft. 

Hence, due to the boundedness of the sequence {p^} C L 2r >^(0, T; W 1 ' i (Cl)) given by 
(16) and (17), along with the almost everywhere convergence of the sequence, we may 
apply the Lebesgue's dominated convergence in order to obtain strong convergence in 
the asserted space. □ 

With the preceding results, we may now pass to the limit r — > in (12). 

Proof of Theorem 1. Let us first establish a weak convergence for {F( r V r )} and 
{$( T )}. Due to estimate (18) we obtain \ € L 2 (0, T; L 2 (Q)) such that 

F (r) p (T) m L 2 (0,T;L 2 (ft)). 

for some subsequence (not relabeled). Due to the almost everywhere convergence for 
p( r ) given in Lemma 11, we may define F such that Fp = \ almost everywhere, where 
we set F(t, x) = when p(t, x) = 0. As for {$( T )}, we obtain from Proposition 3 the a 
priori estimate ||$||2,2 < c ll n ll2, which directly gives us the uniform bound for {<I>( r )}. 
Therefore, there exists some subsequence and $ such that 

$("0 _a $ i nJ L 1+<5 (0,T;I/ ). 

Now, by multiplying (12) with the appropriate test functions, integrating over time 
and space, integrating by parts, and passing to the limit r — > 0, we finally obtain 



(d t n, Vl ) = ((Fp)2Vp, V<pi) + ((Fp)p, A<pi) V<pi e i 9 (0, T; V), 
(Fp,tp 2 ) = (-e 2 Ap + 2p log p-ap*,^) V<^ 2 e L 2 (0, T; L 2 (ft)), 
($, <p 3 > = <«, Vs) V<p 3 e fl^n), 
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with q = {1 + 8)/ 8, which completes the proof. □ 

Appendix A. Without loss of generality, = 0. Otherwise, we may simply 
make the shift vl = u — € Wo and proceed with the proof for u'. From Poincare's 
inequality 

||w-|4t / udx\\ p < Cp||Vu|| p for any u e W 1 '?, 
Jn 

we deduce that the norms |j • \\w arL d || • 1^ _a_ are equivalent. 

'd—1 

We know from standard elliptic theory that a unique weak solution u e Wq of 
(3) exists when / £ L 2 (fl) due to the Lax-Milgram theorem, the Poincare inequality, 
and the continuous embedding i? 1 (fi) =-> Wo- Since L 2 (fl) =-> Llog-L(f2) is dense, it 
is sufficient to show, due to the BLT theorem, that the solution operator S: f \-+ u is 
bounded with respect to the norms || • ||i,iogi,(n) an d II 1 \\w () respectively. 

For (j){s) := sign(s) log(l + \s\) we set 

$(«):= (l + |*|)(log(l + |a|)-l), and 9{ S ) := 2{1 + \s\)^ 2 . 
We define 4 f (r) as the conjugate convex function of $(s), i.e., 

*(r) := sup (rs - $(s)) , 

s 

where the supremum is attained if and only if r = 4>{s) [18]. We first observe that 
|u| € W '' 3=I (Sl) for any u e W 1 ' 3 ^ 1 ^) [17]. Moreover, we have the estimate 

IIVMH^ < ||i + | M |||^||v^)|| 2 < 5 -\\u\\ A + ^||V0(«)||2 + 5 -\n\^i d , 

for any 8 > 0. Clearly, this follows directly from Holder's and Young's inequality. Due 
to the Sobolev embedding W 1 '^ 1 (Cl) <—> L^ r5 (f2) as well as the norm equivalence 
II • 1 1 w ~ II ' Hi d ■ we ma y choose an appropriate 8 > to further obtain 

(19) l|V|u|||_d_ < c\\\79{u)\\l + l-\Sl\( d - 2 )/ d for any we W , 

d -! 2c 

for some constant c > 0, depending only on d and |0|. Now, by testing equation (3) 
with <fr(u), we obtain 

a [ \V8(u)\ 2 dx < [ $(fc/)dx + f V(l<f>(u))dx, 
Jn Jn Jn * 

which holds for any k > 0. Since $ and ^ are convex with $(0) = \£(0) = 0, we have 
$(fcs) < fc$(s) and V (k' 1 <fi(s)) < fc- 1 *(0(s)) = fc-^s^s) - $(s)). Together with 
(19) we get 

Nko<- / $(/)d.T+4 r / ((l + | U |)-l0g(l + | U |))dx+^|0|(' i - 2 )/ d . 

Using the Sobolev embedding W 1,a=T (fi) <-» L 1 ^) and choosing k > appropriately 
finally leads to the estimate 

\\Sf\\w = Nk < c(||/|U logi(n) + M(\Q\)) for all / e L 2 (Q), 
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where c > is a constant depending only on d, |0| and a, and M:K + — > R + is a 
monotonically increasing function of |f2|. Consequently, we may extend the solution 
operator S to a solution operator S: L log L(Q) — > Wq, which establishes the existence 
of solutions for (3). Uniqueness follows by standard arguments involving linearity and 
the maximum principle, and therefore concludes the proof. □ 

Appendix B. Let ^/n E H 2 {Vt) and n E L°°(f2) be uniformly bounded away 
from zero. Then for any Si, <5 2 > and some a E (0, 1) for d = 2 or a E (1/3, 1) for 
d = 3, there exists a constant c > such that 



/ n 2 dx < Si / n 
A2 ./n 



dx + S2 n 



Jo 



dx + c((5i, S 2 ,a, \\n\\i), 



where the constant c depends also on the constants from Proposition 2. 
Indeed, if we set p = y/n, then from the Proposition 2, we obtain 



IHII = IHl4<ci|H|f j2 |H| 2 * - 

On the other hand, by setting z — ^fn, we obtain 



cil|p||? 2 IMIi V - 



n\\o = \\z\ 



<C2\\Z\\1 A 



= C2||2||i, 4 ||«lll 



8 d 
32 



From the two inequalities above, we have 



M| 2 < c ?||HlLlMli V , ll«|| 2 <^||z|ltll«lli"- 



Now let a E (0, 1) for d = 2, a E (1/3, 1) for d = 3, and 



IMlS- a <c? (1 - a) l 



(l-°)d (l-c)(4-d) 

' Inll ' 



12 — c 2 

Putting them together yields, 



IIHI1,2 II "111 
i|l+a 4(l+a)„ n 11 ^,, 11 



16 — (3 — a)d (l + a)d (l-a)d 



2 ll'lll Iklll^ 2 IHIl,2 

< ^llkll 1,4 + <y Hli, 2 + c(Si,5 2 ,a, \\n\\i) 
= Si\\Vz\\i + S 2 \\Vp\\ 2 . + c(Si,S2,a,\\n\\i) 

= Si [ n- 1 \V^i\ 4 dx + S 2 [ |Vv^| 2 dar + c(Ji,<52,a,||n||i), 

which is the required inequality. □ 

Appendix C. In order to show the assertion, we employ the method introduced 
in [15] (c.f. [14]). We introduce the functions 

e = ^, v= 1 7j ^, (,+M)e 2 = 4vV(vp) 2 , 

p d p p A 



where A : (b) 2 = £i,j=i a^b^ for A = (a y ) E R dxd , b E R d , and g > by 

2 

= d V 2 



V 2 p 



rf 



d- 1 



2 , 2 
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which exists due to [15, Lemma 3]. Now set the functionals 



J(P) 



V 2 p 


\^ 


Vp 


2 

-6 


Vp 


4" 


p 


p 


p 




P 





dx 



= / 9 
In 



d V 2 + 



d-1' 



p 2 + 9 2 )+d V e-sz 4 



dx, 



dr, 2 



d-V 



g dx, 



UP)= I |V 2 p| 2 dx= / 
Jn Jn 

and the dummy integral expression 

l(p)= f dW(p- 1 \Vp\ 2 Vp)dx= f p 2 ((d + 2) V e + 2pe-t 4 )dx, 
Jn Jn 

which is zero due to the boundary condition d„y/n = on T. The goal is to find 
constant Co, ci such that J — co/C = J — c JC + c{I > 0, and in such a way that c is 
chosen to be as large as possible. Summing the functionals up as described, we obtain 



(20) (J 



c K + Cl l){p)= [ p 2 [(l-c )dr 1 2 + (d + c 1 (d + 2))r ] e + Q(p,^g)]dx, 
Jn 



where Q is the polynomial in p, £ and g given by 

Q(p, e, Q) = (1 - co)^^ 2 + 2c lM e 2 - (ci + 6)£ 4 + (1 - c )g 2 . 

By choosing c\ = —d/ (d + 2), the second term on the right hand side of (20) vanishes. 
We write Q(p, £, g) = bip 2 + 2b 2 p^ 2 + &3^ 4 + big 2 , where 



h = 



(1 -co)d 
d-1 ' 



bo = 



d+2' 



d+2 



-6, bi = (1 - c ), 



and demand that Q > for any given p, £, g. If c < 1, then 64 > 0. Now we choose 
Co < 1 in such a way that the remaining terms in Q is nonnegative as well. This is 
the case if b\ > and b\b^ — b\ > 0. The second condition is equivalent to 



(1 - c )(d + 2)(1 - {d + 2)5 /d) - (d - 1) > 0, 



and therefore 



c < 1 - 



d-1 



(d + 2)(l - (d + 2)S/d) 



< 1 



d-1 



d + 2 d + 2' 



for S chosen sufficiently small. Choosing such a c € (0,3/(c? + 2)), and using again 
the inequality ||V 2 /9|| 2 < ||Ap|| 2 , we obtain 



{J-colQ(p)> [ p 2 {l-c»)dr, 2 dx= 1 —^ ( |Ap| 2 dx>i 
Jn « Jn 



- c 



IC(P), 



which yields further, J{p) > -fK.(p), with 7 = (1 + (d — l)c )/d. □ 
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